Comparing cavity and ordinary laser cooling within the Lamb-Dicke regime 
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Cavity-mediated cooling lias the potential to become one of the most efficient techniques to cool 
molecular species down to very low temperatures. In this paper we analyse cavity cooling with single- 
laser driving for relatively large cavity decay rates k and relatively large phonon frequencies v. It 
is shown that cavity cooling and ordinary laser cooling are essentially the same within the validity 
range of the Lamb-Dicke approximation. This is done by deriving a closed set of rate equations and 
calculating the corresponding stationary state phonon number and cooling rate. For example, when 
V is either much larger or much smaller than k, the minimum stationary state phonon number scales 
as /16v^ (strong confinement regime) and as k/Av (weak confinement regime), respectively. 

PACS numbers: 03.67.-a, 42.50.Lc 



I. INTRODUCTION 

One motivation for cooling particles by coupling them 
to an optical cavity is the hope to cool many of them 
simultaneously and very efficiently to very low temper- 
atures [J, S]. Beyond this, cavity cooling might open 
the possibility to transfer not only atoms and ions but 
also molecules to very low phonon numbers 0-0] ■ In 
general, molecules do not possess the closed transitions 
required for laser cooling [9l-[T2|. Using for example 
ordinary laser cooling, the spontaneous emission of pho- 
tons would rapidly populate states where the particles 
no longer see the cooling laser. Cavity cooling over- 
comes this problem by depopulating the excited elec- 
tronic states of molecules in a two-step process (l3j . 
First, the population in the excited state is transferred 
into the cavity mode. From there it leaks into the envi- 
ronment via spontaneous emission. 

Indications that cavity couplings might enhance the 
confinement of trapped particles have been observed as 
early as 1997 in Grangier's group in Paris (T^]. System- 
atic experimental studies of cavity cooling have subse- 
quently been reported by the group of Rempe p^l - [l8| . 
Vuletic [l9l - [2l| . and others [22l.l23j. Recent atom-cavity 
experiments access an even wider range of experimen- 
tal parameters by using optical ring cavities [U [2^ and 
by combining optical cavities with atom chip technology 
261 [27| , atomic conveyer belts [1^, H^, and ion traps 
30|. However, these experiments are still limited to the 
cooling of atoms and ions. Moreover, when it comes to 
the cooling of large numbers of particles, evaporative and 
sympathetic cooling still seems a much more efficient ap- 
proach [ni-iii. 

Theoretical work on cavity cooling of strongly confined 
particles has been pioneered by ZoUer's group in 1993. In 
Refs. H^H^, they calculate the cooling rate and the final 
temperature of a single two-level particle trapped inside 
an optical resonator with the help of a quantum optical 
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FIG. 1: Experimental setup of a single two-level particle in- 
side an optical cavity with coupling constant g and sponta- 
neous decay rates k and F. The motion of the particle orthog- 
onal to the cavity axis is strongly confined by an externally 
applied harmonic trapping potential with phonon frequency 
V. The cooling of this vibrational mode can be done with the 
help of the cooling laser with Rabi frequency fi. 



master equation within the Lamb-Dicke regime. A simi- 
lar master equation approach to cavity cooling has been 
used in Refs. 0, [13, [11]. Cavity cooling of free parti- 
cles was first discussed in Refs. [H, H^]. Later Ritsch 
and collaborators d [H EH-Iii^uletic et al. 0, E^, 
Murr et al. [46l - l48j , and others |49| developed semiclassi- 
cal theories to model cavity cooling pro cesses, including 
the cooling of polarisable particles [50l - l52| . Moreover, 
Xuereb et al. [53| introduced a simple input-output for- 
malism which can in principle be applied to a variety of 
cooling scenarios. 

This paper summarises and generalises our current 
knowledge about cavity cooling within the so-called 
Lamb-Dicke regime [l, Issl - lssl . lisj and compares cavity 
cooling with ordinary laser cooling. Moreover, it intro- 
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FIG. 2: Level configuration showing a single particle with 
ground |0) and excited state |1). Here lJc and ojl are the 
frequency of the cavity and of the cooling laser. The cor- 
responding detunings with respect to frequency of the 0-1 
transition, ojo, are A and A + <5. 



duces a methodology which paves the way to an analy- 
sis of cavity cooling beyond the Lamb-Dicke approxima- 
tion [s^ • The setup considered in this paper is shown in 
Figure [T] and consists of a single two-level particle with 
ground state |0) and excited state |1) trapped inside an 
optical cavity. In this paper, we assume confinement of 
the motion of the particle in the direction of the cooling 
laser which enters the setup orthogonal to the cavity axis. 
The aim of the described cooling process is to minimise 
the mean phonon number in this vibrational mode. 

The energy levels considered in this paper arc shown 
in Figured] In the following, we denote the detuning of 
the cavity and of the laser with respect to the 0-1 tran- 
sition of the particle by A and A -I- (5, respectively. Here 
we are especially interested in the case, where A is much 
larger than all other system parameters. In this case, the 
particle remains predominantly in its ground state and 
its electronic states can be adiabatically eliminated from 
the time evolution of the system. As we shall see be- 
low, spontaneous emission from the excited atomic state 
remains negligible throughout the cooling process, when 
the experiment operates well within the so-called strong 
coupling regime, i.e. when the cavity coupling constant 
is much larger than the product of the spontaneous decay 
rates k and F. Outside, the strong coupling regime, the 
cooling process illustrated in Figure [1] becomes a mixture 
of cavity and ordinary laser cooling ^45i] . 

In the following, we use the quantur n op tical master 
equation for atom-cavity systems [35l . l36j to derive a 
closed set of rate equations. These equations are lin- 
ear differential equations which describe the time evolu- 
tion of expectation values. Only two of the variables 
in these cooling equations arc populations: the mean 
phonon number m and the mean number of photons in 
the cavity n. All other variables are coherences. More 



concretely, we assume that the mean phonon number m 
evolves on a much slower time scale than all the other ex- 
pectation values which are included in the cooling equa- 
tions, as it applies for a very wide range of experimental 
parameters. This allows us to simply reduce the above 
mentioned cooling equations to a single effective cooling 
equation via an adiabatic elimination of all expectation 
values other than m. As a result, we obtain the cool- 
ing rate and the stationary state phonon number as a 
function of the system parameters. 

Different from Refs. j36H38| we do not restrict ourselves 
to the case where the detuning 5 is close to the phonon 
frequency v. Here we give analytical results for the sta- 
tionary state phonon number m^^ of the trap ped particle 
for the general case. Different from Refs. [ll. Il3l. [20l . IZTI - 
[52j . we avoid semiclassical approximations. Different 
from Ref. 0, we avoid the rotating wave approximation 
with respect to the phonon frequency v. This allows us 
to calculate cooling rates and stationary state phonon 
numbers correctly even when the phonon frequency v 
is relatively small. In fact, it can be shown that the 
counter-rotating terms in the particle-phonon interaction 
can heat a system very rapidly [55| . 

The general idea behind cavity cooling is the contin- 
uous conversion of phonons into cavity photons. While 
phonons have no decay rate, photons leak out through 
the cavity mirrors easily. The conversion of phonons 
into photons can hence result in the constant removal of 
phonon energy from the system. The role of the cooling 
laser is to establish an effective coupling between the vi- 
brational states of the particle and the cavity field. Cool- 
ing, i.e. the permanent loss of a phonon, occurs when the 
absorption of a phonon is accompanied by the creation 
of a photon in the cavity field which subsequently leaks 
into the environment. 

As we shall see below, in the validity range of the so- 
called Lamb-Dicke approximation, there are many sim- 
ilarities between cavity and ordinary laser cooling. For 
example, it will be shown that one should choose the 
effective laser detuning 



(1) 



equal to v, when the cavity decay rate k is much smaller 
than the phonon frequency v (strong confinement regime) 
in order to minimise the final number of phonons in 
the system. The corresponding stationary state phonon 
number rw^^^^^ will be found to be approximately given 
by K^/ldv^. However, if k is much larger than v 
(weak confinement regime), one should choose ^off — 
since mf i equals n/Av to a very good approxima- 

o^ff — 2 ^ 

tion. Comparing the explicit analytical expressions for 
mf 1 and rnf ,,, we find that 

(5et£ = 5K Ooff = l" 



(5cff — 1/ 



(2) 



for a wide range of experimental parameters. This means, 
when the trapped particle cannot be so strongly confined 
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that it is not possible to cool to phonon numbers below 
one, it is better to choose 5cS = than choosing (5cff = 
v. since rrf^ , < vvf _,, in this case. 

Comparing the explicit analytical expressions for the 
cooling rates 75^ff=i/ and ^s^a^^n fo'" the two above men- 
tioned cases, i.e. the strong and the weak confinement 
regime, we will find that 

to a very good approximation. This means, choosing 
^cs ~ instead of i5cff = v yields a speed up of the 
cooling process when k > 8;^. As we have seen in the pre- 
vious paragraph, this detuning is also the better choice 
in order to minimise the stationary state phonon num- 
ber in the weak confinement regime. Speeding up the 
cooling process is important, since spontaneous emission 
from the excited state |1) remains thus negligible for a 
much wider range of single particle-cavity cooperativity 
parameters / kT . 

There are five sections in this paper. In Section HIl we 
derive the master equation for the experimental setup 
shown in Figure [TJ The main approximations made in 
this section are the Lamb-Dicke approximation, the ro- 
tating wave approximation with respect to optical tran- 
sitions, and the adiabatic elimination of the electronic 
states of the particle. Section IIIII analyses the cavity 
cooling process. Starting from the master equation, we 
obtain a full set of cooling equations which are then used 
to calculate the stationary state phonon number rrf^ and 
the cooling rate 7 as a function of the experimental pa- 
rameters. Section [IVI compares cavity cooling with ordi- 
nary laser cooling. Finally, we summarise our results in 
Section [V] There we emphasize that going beyond the 
Lamb-Dicke regime might yield very different results. 



II. THEORETICAL MODEL 

Let us start by introducing the theoretical model for 
the description of the experimental setup in Figure [TJ It 
contains a strongly confined particle coupled to an optical 
cavity. The cooling laser drives the electronic transition 
of the particle and enters the setup in a direction orthog- 
onal to the cavity axis. In the following, we assume that 
an external trapping potential confines the particle in the 
laser direction such that its motion becomes quantised. 
The aim of the cavity cooling scheme discussed in this pa- 
per is to minimise the number of phonons in this mode. 
Cooling other vibrational modes of the particle motion 
would require additional cooling lasers. 



A. The Hamiltonian 

The Hamiltonian of the system in Figure [T] can be writ- 
ten as 

H = Hpai + -ffphn + Hcav + Hi, + -ffpai-cav • (4) 

The first three terms are the energy of the electronic 
states of the trapped particle, its quantised vibrational 
mode, and the quantised cavity field mode. Suppose, the 
particle is effectively a two-level system with ground state 
|0) and excited state |1) and the energies fkuo, hv, and 
hhJc are the energy of a single atomic excitation, a single 
phonon, and a single cavity photon, respectively. Then 

iJpar = ftwo (T^a~ , 

i/cav = ftWcC^C, (5) 

where the operators = |0)(1| and (t+ = |1)(0| arc the 
atomic lowering and raising operators, and b and c are 
the phonon and the photon annihilation operators with 
the commutator relations 

[6,6t] = [c,ct] = l. (6) 

The two remaining terms in Eq. i.e. Hi, and _ffpar-cav, 
are the laser Hamiltonian and the Hamiltonian describ- 
ing the interaction between the trapped particle and the 
cavity mode. Let us now have a closer look at these 
terms. 

The role of the cooling laser is to establish a coupling 
between the electronic states |0) and |1) of the trapped 
particle and its quantised motion. In the dipole approx- 
imation, its Hamiltonian can be written as 

Hi, = eD-EL(x,t), (7) 

where e is the charge of an electron, D is the atomic 
dipole moment, and EL(x,t) is the electric field of the 
laser at the position x of the particle relative to its equi- 
librium position R = at time t. More concretely, the 
dipole moment D can be written as 

D = Doia-+H.c. , (8) 

where Dqi is a 3-dimensional complex vector, and 

El(x, t) = Eo e'C'L-^-"^*) + c.c. (9) 

with Eq, kL, and ojl being the amplitude, the wave vec- 
tor, and the frequency of cooling laser. 

As already mentioned above, the cooled motion of the 
trapped particle is its center of mass motion in the laser 
direction. Considering this motion as quantised with the 
phonon annihilation operator b from above, yields 

kL-x = + (10) 
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where the Lamb-Dicke parameter ry is a measure for the 
steepness of the trapping potential [l^. Introducing the 
particle displacement operator D as 



(11) 



and substituting Eqs. (|8|)- (fTT|) into Eq. ([7]), we see that 
the laser Hamiltonian can also be written as 

i?L = e [Doicr- +H.C.] •E*i:>(i?/)e'"L* +H.C. (12) 

The cooling laser indeed couples vibrational and elec- 
tronic states of the particle. 

The interaction Hamiltonian describing the coupling 
between the electronic states of the trapped particle and 
the cavity equals in the usual dipole approximation 



par— cav 



= eD • Ecav(x) , 



(13) 



where Ecav(x) is the observable for the quantised electric 
field inside the resonator. Denoting the coupling con- 
stant between the particle and the cavity field as g, this 
Hamiltonian can be written as 



i?par-cav = hg{a +0-+)c + H.C. 



(14) 



which describes the possible exchange of energy between 
atomic states and the cavity mode with a constant cavity 
coupling constant g. 



B. Interaction picture 

We now change into an interaction picture which allows 
us to take advantage of the fact that the experimental 
parameters f2, 6, v, g, and A are in general much smaller 
than the optical frequencies, i.e. 

S , v , g , A < WL, Wc. (15) 

Choosing 

Hq = hbJi^ (7^ + h {ujc — S) c , (16) 

we find that the interaction Hamiltonian Hi, 

Hi = U^o{t,0){H -Ho)Uo{t,0), (17) 

contains terms which oscillate with frequencies close to 
2a;L, 2u!c, ^^l + oJc- Neglecting these relatively fast oscil- 
lating terms as part of the usual rotating wave approxi- 
mation, Hi becomes 

Hi = ^hnD{b])(7^ + hga^c^ +II.C. 



+h{A + S)(T+(T- +hi^b^b + h6c''c, (18) 



where 



2e 



■ Doi ■ Eq 



(19) 



denotes as usual the laser Rabi frequency. This Hamilto- 
nian is time- independent and will be used later to analyse 
the time evolution of the mean phonon number. 



C. Adiabatic elimination of the electronic states 

To minimise spontaneous emission from the excited 
electronic state of the trapped particle, i.e. to maintain a 
closed cooling cycle, we assume in the following that the 
detuning A is much larger than all other system param- 
eters. 



A:^n,s,v,g,r,K. 



(20) 



Concentrating on this parameter regime allows us to 
eliminate the particle adiabatically from the time evo- 
lution. To do so, we write the state vector of the system 
as 



1 oo 
j—O rn,Ti— 



(21) 



where \j) denotes the atomic state and \m) and \n) are 
phonon and photon number states, respectively. Accord- 
ing to the Schrodinger equation, the time evolution of the 
coefficient Cj'm'n' is given by 



-\Y.Y. Cj.m,,{j'm'n'\Hi\jmn) . (22) 



j — O 771,71 — 



Given condition f[2Q|l. the Cj'm'n' with j' = 1 evolve on a 
much faster time scale than the coefficients with j' — 0. 
Setting their derivatives equal to zero, we find 



Clm'71^ 2^ ^ ^ ^077171 

771,71 — 

x(mV| [2g c + n D'' {u-j)] \mn) (23) 



which holds up to first order in 1/A. Substituting this 
result into the differential equations for com'n', we obtain 
the effective interaction Hamiltonian 

Hi = hgcsD{iri)c + R.c. + hiyb''b + hScsc''c (24) 

with the (real) cavity coupling constant goff given by 

(25) 



Soft 



gn 

2A 



and the effective detuning 5cS as defined in Eq. The 
interaction Hamiltonian Hi in Eq. holds up to first 
order in 1/A. It no longer contains any atomic operators. 
Instead. Eq. describes a direct interplay between 

phonons and cavity photons. 



D. Lamb-Dicke approximation 

Suppose the particle has already been cooled enough 
to ensure that it remains in the vicinity of its equilib- 
rium position R = 0. More concretely, we assume in 
the following that the displacement x of the particle is 
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small compared to the wavelength of the cooling laser. 
Then Icl • x in Eq. (jlOp is much smaller than one, and 
the so-called Lamb-Dicke approximation with 

V < 1 (26) 

can be applied. This means, Eq. simplifies to 

D{ir]) = 1-iriib + b''). (27) 



Substituting this into Eq. p4|. we finally obtain the in- 
teraction Hamiltonian 

Hi = hgcs c ~ ihrjgcff (b + b^)c + }i.c. + hvb^b 

+Mcffctc (28) 

which contains cavity interactions, the phonon-photon 
interaction, the phonon energy term, and a level shift. 



E. Master equation 

After the adiabatic elimination of the electronic states 
of the particle, the only relevant decay channel in the 
system is the leakage of photons through the cavity mir- 
rors. To take this into account, we describe the cooling 
process in the following by the master equation in the 
usual Lindblad form 



P 



UHi,p] + 1-k(2 cpc'^ - c)cp - pc)c) (29) 
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with Hi as in Eq. ^ 



III. ANALYSIS OF THE COOLING PROCESS 

In this section, we calculate the stationary state 
phonon number rrf^ and the efi^ective cooling rate 7 as 
a function of the experimental parameters 77, gos, k, z^, 
and 5eS- This is done by using the above master equation 
to obtain a closed set of cooling equations which predict 
the time evolution of the mean phonon number m. In the 
following, we are especially interested in the case where 
the cavity decay rate k and the phonon frequency v are 
both relatively large. More concretely, we assume in the 
following that 



(30) 



As we shall see below, the mean phonon number m 
evolves in this case on a much slower time scale than 
all other expectation values which are included in the 
cooling equations. The latter can hence be eliminated 
adiabatically from the time evolution of the system. 



A. Cooling equations 

Since the time derivative of the expectation value of 
an operator A equals 



the master equations in Eq. (|29|) implies 

(i) = {[A,Hi]) + ]^k{2c'^Ac- Ac^c-c^cA) .{il) 

Using the commutator relations in Eq. ^ and the inter- 
action Hamiltonian in Eq. (j28p and applying this equa- 
tion to the mean phonon number m = (6^6) and the mean 
photon number n = (c^c), and the coherences 

fc, = i(5-&t), fcy=i(C-ct), 

fcu = (6 + 6^), fcw = (c + ct), 

fci = {{b + 6t)(c + ct)} , fc2 - i((& + b^){c - ct)) , 

h = i((6 - b^){c + ct)) , h = {{b - b^){c - ct)) , 

h = {b^ + 6^2) , ks = i(&2 - 6t2) (33) 
we obtain a closed set of differential equations. These are 
ftx = -^rjgcs ky + vku, 



ky = 2gcs + Scs k^ - -Kky , 



ku — V k-y^ . 



1 



few = 2?75off' fcu - ^csky - -K k. 



(34) 



and 

n = gcS ky + rjgcs fci - k n , 

fci = 27/goff (fcy 4- 2m + 1) - J/ A:3 - S^s k2 - ^nki , 

k2 = 2gcfi ku + 1^ k^ + 6cs ki - k2 , 

fcs -2ry5off (ke ~ kg) + ly ki + S^fi k^ ~ ^Kk^ , 

ki = -2gcff fcx - 2775cff (^5 - 2n - 1) - A:2 - 4ff ^3 
--Kki, 

k5 = -2goff ky + 2r]gcs ki - 2(5cff kg - nk^, 

h = 25off few + 2?7goff k2 + 2ScS k^- nk^, 

kr = -2?7goff ^4 ~2iyks, 

ks = -2r]gcs k2 + 2iy kr , (35) 



while 



(36) 



{A) = Tr{Ap), 



(31) 



Notice that these differential equations, the cooling equa- 
tions, have been derived without further approximations. 



B. Stationary state phonon number 

Calculating the stationary state phonon number rrf^ 
can now be done by setting the right hand side of the 



6 



above cooling equations equal to zero. Doing so, Eq. ((341 
yields for example 



0, 



with the cubic frequency defined as 
Moreover, we obtain the stationary state values 



(37) 



(38) 



kf 



f^ss _ 



V9c«k{k'^ + 4(5gff) 



cff 



vgcff 



ss 



4 - 



4<S, 



0, 



off 



) , ri^gl^{^' 



4<Se'ff) 



7,SS 



32ffgffKi^^(5cff 



4??^3cffK 



(39) 



and 



m 



, K«2 + 4J2^)(^ 



16i/(5, 



cff 



452 



off 7 



2(5, 



'off J 



, 16^^5o^ff-^ 



45effM3 ■ ■ 

These equations can be checked easily by substituting 
them back into Eqs. 

In this paper we are especially interested in the pa- 
rameter regime of a tightly confined particle inside a rel- 
atively leaky optical cavity described by Eq. (PH)) . This 
parameter regime is consistent with the Lamb-Dicke ap- 
proximation in Eq. p6p . Calculating m^'^ up to second or- 
der in rj^ correctly would require to go beyond the Lamb- 
Dicke approximation and to take terms proportional 77^ 
in the system Hamiltonian into account. The above ex- 
pression for the stationary state phonon number hence 
applies only up to first order in rj. Taking this into ac- 
count, Eq. (jini) simplifies to 



m 



16i/(5off 



(41) 



0.08 



0.06 - 




0.04 



0.02 - 



FIG. 3: Logarithmic contour plot of the stationary state 
phonon number Tif in Eq. (|40|l as a function of Soft and u 
for 77 = 0.1 and goff = 0.0001 k. The result is in very good 
agreement with the simpler expression in Eq. (|4ip . 



Since Eq. (pO| does not pose a condition on the size of the 
effective detuning Joff, this parameter can now be used 
to minimise the stationary state phonon number rrf^ in 
Eq. dHI). 

Calculating the derivative of m^^ with respect to ScS 
we find that the optimal choice for (5off is 



OofT 



5^ 



4l^2 



(42) 



In the parameter regime which is the most interesting 
from an experimental point of view, i.e. in the case of 
relatively small phonon frequencies v (weak confinement 
regime), the effective detuning in Eq. (|^^ becomes 



ScS 



(43) 



Substituting this detuning into Eq. (|41l) . we obtain the 
stationary state phonon number 



TO. 



4i^ 



(44) 



which is in good agreement with Figure [31 In the other 
extreme case, i.e. when the phonon frequency v is much 
larger than the cavity decay rate k (strong confinement 
regime), the effective detuning in Eq. (|42|) simplifies to 



^off = f 

which corresponds to 

ss 



(45) 



(46) 
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100 




100 



FIG. 4: Contour plot of the stationary state phonon number 
in Eq. (|40|) as a function of S^ff and for the same rj and 
geff as in Figure [3] but for much larger phonon frequencies i/. 
The result is again in very good agreement with Eq. (|4ip . 



This result is confirmed by FigureUwhich considers much 
larger phonon frequencies v than Figure [S] 

Comparing the two choices of effective detunings 6cs 
in Eqs. P5|) and we find that cooling to very low 

temperatures is always limited by the relative size of the 
phonon frequency v with respect to the cavity decay rate 
K. However, comparing Eq. (|4T1) and Eq. (|44|) . we see 



that the stationary state phonon number for Sctf = is 
the square root of the stationary state phonon number 
for 6cff = v, as it has already been pointed out in Eq. ([2|l. 
This means, choosing the detuning Scs different from the 
phonon frequency as it has been suggested in Refs. (36l - 
Issl . |45| . yields a significant enhancement of the cavity 
cooling process exactly when it is most needed, namely 
when m^^ is larger than one. 



for a = X, y,u, w and i = 1,...,8, while there can be 
any mean number of phonons m in the vibrational mode 
of the particle. This initial condition is consistent with 
the particle being trapped which means that it is located 
around the centre of a trap and that it has no initial 
momentum away from its equilibrium position. The first 
of these two statements implies fcu(O) = and the second 
one implies fcx(O) = 0. 



D. Cooling dynamics 

To calculate the effective cooling rate 7, we notice that 
only one of the variables in the above cooling equations, 
namely the mean phonon number m, evolves on a rela- 
tively slow time scale. All other variables, i.e. the mean 
photon number n and the coherences, evolve on the fast 
time scale given by v and n. In the parameter regime 
of Eq. (PH)) . these can hence be eliminated adiabatically 
from the time evolution of the system, leaving us only 
with a single effective cooling equation. Doing so and 
setting the time derivatives in Eqs. (|34p and ([55)1 equal 
to zero and assuming that we are at the beginning of the 
cooling process where m ^ 1, we find 



^4 = — 



(48) 



This equation holds up first order in Tyt/cff and is consis- 
tent with the Lamb-Dickc approximation introduced in 
Section IUdI 

Eq. p6l) shows that the photon-phonon coherence is 
essentially the cooling rate of the trapped particle. Sub- 
stituting Eq. (pS]) into Eq. ([55]) . we obtain the final cool- 
ing equation 



m = — 7TO 
with the cooling rate 7 given by 



1 



Solving this equation yields 



(49) 



(50) 



Initial state 



m{t) 



-'ft 



to(0) . 



(51) 



To determine the initial state of the experimental setup 
in Figure m we assume that the cooling laser is turned on 
at t = 0. Moreover, we assume that the particle does not 
experience any other cooling processes and is at t = in 
an equilibrium state of the master equation ([29| but with 

= 0. Taking this into account and setting gcff and the 
right hand side of the cooling equations (|34)) -(p6 | equal 
to zero, we find that this corresponds to a state with all 
coherences and the cavity photon number being equal to 
zero, i.e. 

n(0) = fca(0) = fc»(0) = (47) 



Since the rate 7 is always positive, this equations de- 
scribes an exponential reduction of the initial phonon 
number 77i(0), i.e. cooling. Figures [5] and |6] compare the 
exponential cooling process with the rate 7 in Eq. (|50p 
(dashed lines) with an exact numerical solution of the full 
set of cooling equations (solid lines) and shows that 7 is 
indeed a very good approximation for the cooling rate as 
long as the actual phonon number m is much larger than 
one. The exponential reduction of m only slows down 
when m approaches its stationary state value. We also 
see that the speed of the cooling process increases, as v/k 
increases. 
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FIG. 5: Logarithmic plot of the time evolution of the mean 
phonon number m for rj = 0.1, QbS = 0.0005 k, and 5c« ~ 
i K for two different phonon frequencies v. The solid lines 
have been obtained from a numerical solution of the cooling 
equations (|34|| - (|36p for the initial conditions in Eq. (|47|) and 
m(0) = 2500. The dashed lines assume an exponential cooling 
process with the rate 7 in Eq. (|50|l . Both solutions coincide 
very well when m is far away from its stationary state value. 



FIG. 6: Logarithmic plot of the time evolution of the mean 
phonon number m for the same experimental parameters as 
in Figure [S] but with SeS = v- Again, the solid lines have been 
obtained from a numerical solution of the cooling equations 
(l3ll)-((36ll for the initial conditions in Eq. ((47| and m(0) = 
2500. The dashed lines assume an exponential cooling process 
with the cooling rate 7 in Eq. (|50|) . Although k and v remain 
the same, we now observe slower cooling processes with higher 
stationary state phonon numbers nf^. 



As in the previous section, let us now have a closer look 
at the two extreme cases, where the phonon frequency 
y is either much smaller or much larger than the cavity 
decay rate k. In the first case (weak confinement regime), 
one should choose (5cff = \n (c.f. Eq. (031)) in order to 
minimise the stationary state phonon number. Doing so, 
Eq. ([50)1 simplifies to 



k4 + 41^4 



(52) 



In the latter case (strong confinement regime), one should 
choose (Soff = v (c.f. Eq. (1451) ) in order to minimise the 
stationary state phonon number. This choice of im- 
plies 



(53) 



Both cooling rates scale as ij^g^g and are hence in general 
relatively small (c.f. Eq. ([30])). As already pointed out in 
Eq. (|3]), comparing the analytical expressions in Eqs. (|52p 
and (|53|) we find that 75^^^^ is about k/Sv times larger 
than 75^j.,.=iy. This is confirmed by Figures [S] and H] which 
show that choosing Scs = in the weak confinement 
regime not only leads to a lower stationary state phonon 
number but also to a significant speed up of the cooling 
process. Large cooling rates arc important when the pur- 
pose of using a cavity is to avoid spontaneous emission 
from the particle. 



E. Avoiding spontaneous emission from the 
particle 

The analysis of the cooling process presented in this 
paper only applies, when the population in the excited 
atomic state |1) remains negligible. Otherwise, the adi- 
abatic elimination in Section III CI does not hold and the 
cooling process becomes a mixture of cavity and ordi- 
nary laser cooling (45j . Avoiding spontaneous emission 
from the excited electronic state |1} is especially impor- 
tant when it comes to the cooling of molecules, where it 
could result in the population of states, where the parti- 
cle no longer experiences the cooling laser. We therefore 
conclude this section with an estimation of the parame- 
ter regime, where spontaneous emission from the particle 
remains highly unlikely. 

In the Lamb-Dicke limit and the parameter regime 
given by Eq. (pil)) . Eq. (P5|) shows that the population 
in |1) scales essentially as f2^/A^. We therefore assume 
in the following that 



7 > 



4A2 



(54) 



i.e. that the cooling rate is much larger than the proba- 
bility density for the spontaneous emission of a photon 
from the particle. For SeS = \k. and 7 as in Eq. ([5^ 



and when taking the definition of goff in Eq. ([25]) 
account, we see that this condition applies when 



into 



(55) 
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Since 77 ^ 1, the right hand side of this equation is in 
general much larger than one. This means, spontaneous 
emission from the particle is only negligible, when the 
cavity is operated in the so-called strong coupling regime. 
If the cavity decay rate n is much smaller than Av, one 
should choose ScS = v and the cooling rate simplifies to 
the expression in Eq. ([55)1 . In this case, condition ([M)) 
simplifies to 



> 



g4jy2jy2 



(56) 



Since 77 ^ 1, we find again that the cavity needs to be 
operated within the strong coupling regime. 



IV. COMPARISON WITH ORDINARY LASER 
COOLING 

In the parameter regime of a tightly confined parti- 
cle inside a relatively leaky optical cavity described by 
Eq. pop , the cavity cooling scenario in Fig ure[T] has many 
similarities with ordinary laser cooling [9l4l2l|. The rea- 
son is that the atomic O-I transition and the cavity are so 
strongly detuned that the electronic states of the trapped 
particle can be adiabatically eliminated from the system 
dynamics (c.f. Section lll C[) . The remaining master equa- 
tion (c.f. Eq. (123)) with the interaction Hamiltonian Hi 
in Eq. (p8)) is almost the same as in laser cooling. One 
only needs to replace the cavity annihilation operator c 
by the atomic lowering operator |0)(1|, the cavity decay 
rate k by the atomic decay rate F, and the effective cou- 
pling constant geS by the cooling laser Rabi frequency $7, 
and so on. A crucial advantage of cavity cooling com- 
pared to ordinary laser cooling could be the lack of recoil 
heating. However, this difference does not matter as long 
as there is always only a very small number of photons 
inside the cavity. 



A. Comparison of cavity cooling with ordinary 
laser cooling 

As in laser cooling, we found that it is useful to dis- 
tinguish between two different parameter regimes: the 
strong confinement regime and the weak confinement 
regime. In the strong confinement regime, where the rel- 
evant spontaneous decay rate (k or F) is much smaller 
than the phonon frequency one should choose the rel- 
evant laser detuning equal to v (c.f. Eq. (|45p ) in order 
to minimise the stationary state phonon number rrf^. In 
laser cooling, this case is known as laser sideband cool- 
ing. The stationary state phonon number rrf^ for cavity 
cooling in the strong confinement regime is essentially 
given by k^/I&v"^ (c.f. E g. (HTI) ') while scaling as V'^/v'^ 
in laser sideband cooling ^M,- This means, it is possible 
to cool to phonon numbers well below one (ground state 
cooling), although realising k <^ v in cavity cooling is 
experimentally very demanding. 



In the weak confinement regime, where the relevant 
spontaneous decay rate (k or F) is much larger than the 
phonon frequency v, one should choose the relevant laser 
detuning equal to (c.f. Eq. (|13])) or ^F, respectively, 
in order to minimise the stationary state phonon number 
m^^ . In this case, the stationary state phonon number 
rrf^ equals k/Av (c.f. Eq. (jH])) in cavity cooling while 
it scales as T/v in laser cooling [l^]. This is exactly 
what one would expect when the described cavity cool- 
ing scheme and ordinary laser cooling are more or less 
equivalent. 



B. Comparing ScH — v with 5cff = \k va. cavity 
cooling 



Previous papers (see e.g. Refs. [36l - l38l l45l|) mainly fo- 
cus their analysis on cavity cooling in the strong con- 
finement regime, where one should choose (5eff = v and 
where it is in principle possible to cool the trapped par- 
ticle to phonon numbers well below one. The purpose of 
this paper is to point out that there are three distinct 
advantages in choosing (Sofi differently, i.e. close to \k, 



(c.f. Eq. (|42|) ). when it is experimentally not possible to 
enter the strong confinement regime: 

1. A reduction of the stationary state phonon number. 
As already pointed out in Eq. S, nf^ 1 equals 

the square root of the stationary state phonon num- 
ber Ttf^^^^^ and is hence significantly smaller than 
'^^i,ii=v ^ wide range of experimental parame- 
ters. This result is confirmed by Figure [7] which 
shows rn^^ jmf^^^^^ for relatively small phonon fre- 
quencies V and a wide range of effective detun- 
ings (5off. In order to minimise the stationary 
state phonon number, one should choose Jeff as in 
Eq. (gg). 

2. An increase of the cooling rate. Calculating the 
ratio 75„,.,.= i„/75„ff=zy using Eqs. JSH) and ([32]) for 
fixed values of r\, gos, k, and v, we find 



2^ 



4k4 



I6j/4 



(57) 



which scales approximately as k/8v, as already 
pointed out in Eq. ([3]). This means, choosing das 
close to i/i results in a significant speedup of the 
cooling process. This result is confirmed by Figure 
[8] which shows 7/75^^=^ for the same parameters 
as in Figure [71 

3. Minimising spontaneous emission from the excited 
electronic state |I) of the trapped particle. This 
is important, when it comes for example to the 
cooling of molecules, where such an emission might 
populate states, where the particle no longer expe- 
riences the cooling laser. As pointed out in the pre- 
vious paragraph, the cooling rate 7 is much higher 
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FIG. 7: Contour plot of rrf^ /raY^^^=^ as a function of 5cS and 
V. This plot has been obtained using Eq. (|40p for 77 = 0.1 
and gcff = 0.0001 k and shows that choosing the detuning S^s 
comparable to k leads to much lower stationary state phonon 
numbers than when v <^ n. 



FIG. 8: Contour plot of 'y/^s^n^v as a function of ScS and v 
for the same parameters as in Figure [T] This plot has been 
obtained using Eq. (|50|l and shows that that choosing the 
detuning Sch comparable to k leads to a significant speed up 
the cooling process when u <^ k. 



when ^off is close to 2'*- As a consequence, the re- 
strictions which need to be posed on the minimum 
size of the single particle cooperativity parameter 
I kT are therefore much weaker in this case. The 
reduction of the cooling time might moreover help 
to balance unconsidered heating processes which 
are, for example, due to stray fields. 



V. CONCLUSIONS 

This paper analyses cavity cooling of a trapped particle 
with ground state |0} and excited state |1) in the parame- 
ter regime, where the phonon frequency v and the cavity 
decay rate k are both relatively large. As shown in Fig- 
ure [l] we assume that the motion of the particle orthog- 
onal to the cavity axis, i.e. in the direction of the cooling 
laser, is strongly confined. The aim of the cooling process 
is to minimise the number of phonons in this vibrational 
mode. The role of the cooling laser is to establish a cou- 
pling between the phonons and the electronic states of the 
trapped particle. The result is the continuous conversion 
of phonons into cavity photons which subsequently leak 
into the environment. The finite temperature limit of 
the cooling process is due to the counter-rotating terms 
in the system Hamiltonian which simultaneously create 
a phonon and a cavity photon. 

We describe the time evolution of the system using 
the usual master equation. The main approximations 



made in the derivation of the master equation arc the 
Lamb-Dickc approximation, the rotating wave approxi- 
mation with respect to optical transitions, and the adi- 
abatic elimination of the excited atomic state |1) of the 
particle. The second approximation is a standard ap- 
proximation in quantum optics and in general in very 
good agreement with experimental findings. The third 
approximation requires a sufficiently large detuning A of 
the 0-1 transition of the particle with respect to the cav- 
ity field and the cooling laser (c.f. Figured]). Without 
further approximations, the master equation can then be 
used to obtain a closed set of differential equations, the 
cooling equations, which predict the time evolution of the 
mean phonon number. To calculate the stationary state 
phonon number, we set the right hand side of the cool- 
ing equations equal to zero. The cooling rate is obtained 
via an approximate solution of the cooling dynamics with 
the help of an adiabatic elimination which applies for rel- 
atively large k, t^, and m and is in very good agreement 
with numerical simulations. 

Our results show that, in the case of a relatively tightly 
confined particle inside a relatively leaky optical cavity, 
the cavity cooling scenario in Figure [T] is essentially equiv- 
alent to ordinary laser cooling. This applies, since the 
mean number of photons inside the cavity remains neg- 
ligible in the parameter regime considered here. If the 
purpose of the cavity is to avoid photon emission from 
the excited electronic state |1) of the trapped particle, 
as it might be necessary for the cooling of molecules, the 
cavity needs to be operated well within the strong cou- 
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pling regime. Otherwise, the coohng process described 
here simply becomes a mixture of cavity and ordinary 
laser cooling. 

Before concluding this paper we would like to point 
out that the above results which have also been ob- 
tained by other authors are only valid within the 
Lanib-Dicke regime. It is usually believed that the 
Lamb-Dicke approximation, i.e. the expansion of the 
system Hamiltonian as a function of 77, is well justified 
for a sufficiently pre-cooled particle. However, when 
avoiding this approximation, the methods used in this 
paper no longer yield a closed set of rate equations. 
In fact, higher-order terms in 77 which have been 
neglected here could, in principle, evolve the system 
in the long term into a state which is very different 
from the stationary state obtained above [5^. When 



calculating stationary states, a priori assumptions about 
the existence and the yy-dependence of the stationary 
state phonon number should be avoided. Moreover, 
going beyond the Lamb-Dicke regime yields corrections 



to the cooling rate 7 of the order r]^ [37 
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